Mechanical and chemical bonding properties of ground state BeH2 
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The crystal structure, mechanical properties and electronic structure of ground state BeH2 are 
calculated employing the first-principles methods based on the density functional theory. Our cal- 
culated structural parameters at equilibrium volume are well consistent with experimental results. 
Elastic constants, which well obey the mechanical stability criteria, are firstly theoretically acquired. 
The bulk modulus B, Shear modulus G, Young's modulus E and Poisson's ratio v are deduced from 
the elastic constants. The bonding nature in BeH2 is fully interpreted by combining characteristics 
in band structure, density of state, and charge distribution. The ionicity in the Be— H bond is 
mainly featured by charge transfer from Be 2s to H Is atomic orbitals while its covalency is dom- 
inated by the hybridization of If Is and Be 2p states. The valency in BeH2 can be represented as 
ggi.99+ jjo.63- ^ which suggests that significant charge transfer process exists. 

PACS numbers: 68.43.Bc, 68.43.Fg, 68.43.Jk, 73.20.Hb 



I. INTRODUCTION 

There has been a great interest in simple alkali-metal 
and alkaline-earth-metal hydride systems stimulated by 
both the fundamental interests and various applications. 
Such kinds of metal hydrides, e.g., beryllium and magne- 
sium and beryllium-magnesium-based hydrides, which all 
have high weight percentage of hydrogen, are promising 
candidates for hydrogen storage. If they become metal- 
lic when subjected to high pressures, they may be can- 
didates for high-temperature superconductivity. Among 
them, MgH2 has been extensively studied both exper- 
imentally and theoretically. Due to the experimental 
difficulties in material synthesis^, on the other hand, a 
comprehensive understanding of BeH2 is still being in 
progress. It had long been considered that Bell2 was 
a polymer and did not crystallize without ternary ad- 
ditions. Recently, however. Smith et al? have success- 
fully synthesized crystalline Bell2 and first established 
the crystal structure as the body-centered orthorhombic 
by synchrotron-radiation-based powder x-ray diffraction. 
Since this experimental establishment, to date theoretical 
studies on physical properties of crystalline BeH2 are very 
scarce in the literature. Vajeeston et alA and Hantsch et 
al.'^ have studied the structural stability of Bell2 by per- 
forming first-principles total-energy calculations for var- 
ious types of structural variants, among which they con- 
firmed that the experimentally observed Bcll2 modifica- 
tion has the lowest total energy. No further theoretical 
studies have been reported, although more experimental 
works are going on to investigate the pressure dependence 
of physical properties of Bell2^. 

In this paper, we have carried out the first-principles 
calculations of ground-state behavior of Bell2 at its ex- 
perimentally established crystalline phase. Results for 
the electronic and atomic structure, mechanical proper- 
ties, charge density participation, and the nature of Be- 
ll chemical bonding, are systematically presented. Our 



calculated results for elastic constants indicate that the 
experimentally observed orthorhombic phase of BeH2 is 
mechanically stable. We determined that the ionic charge 
of Be and H in Bell2 are represented Be^'^^"*" and H"-^"^", 
which, when compared to Mg^-^^+ and Il"-^~ in MgH2, 
indicates that the charge transfer process in Bell2 is more 
prominent in Mgll2. 



II. CALCULATION METHOD 

The density-functional theory (DPT) total energy 
calculations for the ground state Bell2 were car- 
ried out using the Vienna ab initio simulations 
package (VASP)^ with the projected-augmented-wave 
(PAW) pseudopotentials^ and plane waves. The gener- 
alized gradient approximation (GGA)^ for the exchange- 
correlation potential was employed. 7x7x7 Monkhorst- 
Pack^ k points in the full wedge of the Brillouin zone 
were used. The plane-wave energy cutoff was set 500 
eV, and all atoms were fully relaxed until the Hellmann- 
Feynman forces were less than 0.002 eV/A. To obtain 
optimized lattice vectors and atomic coordinates of bulk 
unit cell, we relaxed it at a series of fixed volumes. The 
bulk modulus B, shear modulus G, Young's modulus E, 
Poisson's ratio u, were gained through computing elastic 
constants. 



III. RESULTS AND DISCUSSION 
A. Atomic structure and mechanical properties 

The ground state of Bell2, which has been testi- 
fied both experimentally^- and theoretically'^, belongs to 
body-centered orthorhombic structure with space group 
Ibam. Its orthorhombic unit cell is composed of twelve 
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TABLE I; Calculated structural parameters, bulk modulus B, shear modulus G, Young's modulus E, Poisson's ratio v and 
energy band gap (Eg) for orthorhombic BeH2. As a comparison, other theoretical work and available experimental data are 
listed. 





Present calculations 


Previous calculations 


Expt. 


a (A) 


9.0049 


8.9823 


9.082 


b{k) 


4.1684 


4.1563 


4.160 


c{A) 


7.6347 


7.6455 


7.707 


V (A3) 


286.6 


285.4 


291.2 


Coordinates 


Bel(4a): 0, 0, 0.25 


0, 0, 0.25 


0, 0, 0.25 




Be2(8j): 0.1684, 0.1199, 


0.1678, 0.1207, 


0.1699, 0.1253, 




Hl(16fc): 0.0881, 0.2216, 0.1516 


0.0720, 0.1842, 0.1361 


0.0895, 0.1949, 0.1515 




H2(8j): 0.3092, 0.2792, 


0.3097, 0.2790, 


0.3055, 0.2823, 


Bel-Hl (A) 


1.328 


1.435* 


1.38(2) 


Be2-Hl (A) 


1.376 




1.41(2) 


Be2-H2 (A) 


1.434 


1.439' 


1.44(2) 


Hl-Bel-Hl (°) 


109.6 




107.7(5) 


Hl-Be2-Hl (°) 


98.3 


98.6* 


112.1(5) 


Hl-Be2-H2 (°) 


117.9 




111.2(5) 


H2-Be2-H2 (°) 


109.2 




109.1(5) 


Bel-Hl-Be2 (°) 


133.6 


126.4* 


130(1) 


Be2-H2-Be2 (°) 


125.4 


120.0' 


127(1) 


B (GPa) 


22.60 


23.79 


14.2±3.0" 


G (Gpa) 


24.50 






E (GPa) 


53.99 






V 


0.102 




0.21" 


Eg (eV) 


5.52 


5.51 





"Experimentally measured values from Ref.— for amorphous BeH2 at ambient pressure. Other 
from Ref.=. 'Previous theoretical values from Ref.* and other from 



experimental values are selected 
Ref.\ 




FIG. 1; (Color online) Unit cell for body-centered orthorhom- 
bic BeH2. Both Be and H have two kinds of atomic occupa- 
tions, which are in the figure presented by blue (Bel), light 
blue (Be2), gray (HI), and white (H2) spheres. 



BeH2 formula units, see Fig. 1. The primitive cell, con- 
taining six formula units, has two Bel (4a, in Wyckoff 
notation), four Be2 {8j), eight HI (16fc), and four H2 
(8j) atoms. Each Bel is surrounded by four HI atoms 
to build the tetrahedral structure and each Be2 connects 
with two HI atoms and two H2 atoms. This can be 



clearly seen from Fig. 1, in which the larger blue spheres 
stand for Bel atoms and the larger light blue Be2 while 
the smaller gray spheres denote HI atoms and the smaller 
white H2. Calculated atomic structural parameters for 
orthorhombic BeH2 are presented in Table I. As a com- 
parison, the previous theoretical results^ji and available 
experimental datai^i^ are also listed. It is clear that our 
calculated results are on the whole in good agreement 
with the experimental data, which thus confirms that the 
present calculation is reliable and accurate. The excep- 
tion is that in analogy with previous theoretical report^, 
our calculated Hl-Be2-Hl bond angle (98.3°) is underes- 
timated by 12% compared with the experimental value 
(112.1°)'^. Although it remains unclear for us to explain 
this discrepancy, we can be sure that it is not caused 
by the choice of the exchange-correlation potential. In 
fact, we have found that the use of local density approxi- 
mation (LDA) even more underestimates the Hl-Bc2-Hl 
bond angle. 

Elastic constants can measure the resistance and me- 
chanical features of crystal to external stress or pressure, 
thus describing the stability of crystals. For the present 
orthorhombic BeH2, there are nine independent elastic 
constants, which are in our first-principles calculations 
obtained by applying small elastic strains. For each time 
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TABLE II: Calculated elastic constant of BeH2 (Ibam). All 



values 


are in 


units of GPa. 












Cn 


C22 


C33 C44 


C55 


Cee 


C12 


C23 


C13 


63.85 


40.33 


89.44 14.69 


34.00 


21.01 


0.63 


0.66 


10.76 



that the strain changes, the coordinates of ions are fully 
relaxed to ensure the energy minimum rule. In our cal- 
culations the strain S is varied in steps of 0.006 from 
i5=-0.036 to 0.036. The calculated elastic constants are 
listed in Table II. Although at present there exist no ex- 
perimental data to compare with our calculated elastic 
constants, based on the good agreement between our cal- 
culations and the experiment for the atomic structural 
parameters, we expect that our calculated elastic con- 
stants can be adopted in the realistic application if the 
mechanical data are needed for BeH2. 

The availability of the elastic constants can be used 
to measure the stability of a specific material. For the 
orthorhombic crystalline structure, the mechanical sta- 
bility criteria are given by^*^ 

Cn > 0, C22 > 0, C33 > 0, C44 > 0, C55 > 0, Cee > 0, 
[Cn + C22 + C33 + 2(Ci2 + Ci3 + C23)] > 0, 
(Cn + C22 - 2C12) > 0, (Cn + C33 - 2C13) > 0, 
(C22 + C33 - 2C23) > 0. (1) 

From Table II we can see that the orthorhombic BeH2 is 
mechanically stable because its elastic constants satisfy 
formula (1). In fact, the calculated amplitudes of C12 
and C23 are negligibly small, and the amplitude of C13 
is prominently smaller than that of either Cn or C33. 
Thus, formula (1) is easily satisfied. 

After obtaining clastic constants, we can calculate 
bulk and shear moduli from the Voigt-Reuss-Hill (VRH) 
approximationiiii^iii^. The Voigt bound s^ ^1^^ on the ef- 
fective bulk modulus By and shear modulus Gy are 

By = (1/9) [Cn + C22 + C33 + 2(Ci2 + Ci3 + C23)] (2) 
and 

Gv = (1/15) [Cn + C22 + C33 + 3(C44 + C55 + Cee) 
-(C12+C13 + C23)]. (3) 

Under Reuss approximation^^, the Reuss bulk modulus 
Bji and Reuss shear modulus Gr are 

Br — A[Cii(C22 + C33 — 2C23) + C22(C33 — 2C13) 

— 2C33C12 -I- Ci2(2C23 — C12) -I- Ci3(2Ci2 — C13) 
+ C23(2Ci3-C23)]~' (4) 

and 

Gr = 15{4[Cn(C22 + C33 + C23) + C22(C33 + C13) 
+ C33C12 — Ci2(C23 + C12) — Ci3(Ci2 + C13) 
- C23(Ci3 -f C23)]/A + 3[(1/C44) + (I/C55) 
+ (l/C66)]}-\ (5) 



where 

A = Ci3(Ci2C23 — C13C22) + C23(Ci2Ci3 

— C23C11) -I- C33(CiiC22 — C^2). (6) 

The bulk modulus B and shear modulus G, based on 
Hill approximation^'^, are arithmetic average of Voigt 
and Reuss elastic modulus, i.e., B—^{Br + By) and 
G—\{Gr -\- Gy). The Young's modulus E and Poisson's 
ratio V for an isotropic material are given byi^ 



35 + G 

and 

3B-2G 
" ^ 2{3B + G)' 

The calculated results for these moduli and Poisson's 
ratio are listed in Table I. Note that we have also cal- 
culated the bulk modulus B by fitting the Murnaghan 
equation of state. The derived bulk modulus is well com- 
parable with that from the above VRH approximation, 
which again indicates that our calculations are consistent 
and reliable. We find that the bulk modulus of BeH2 is 
much smaller than that of MgH2 (previous first- principles 
reports^^'^*' gave 5^50 GPa in the MgH2 system), which 
means that BeH2 is more easily compressed. The rea- 
son is that the rutile structure of MgH2 is more dense 
in atomic packing than the body-centered orthorhombic 
structure of BeH2 . Although at present no experimental 
data for these moduli and Poisson's ratio are attainable 
for pure BeH2, recent Brillouin scattering experiment^ 
on amorphous BeH2 gave 5=14.2±3.0 GPa and u=0.21. 
Considering the fact that the amorphous BeH2 is more 
loose than the pure BeH2, it is understandable that our 
calculated modulus for the pure phase is larger than 
the experimental measurement for the amorphous phase. 
Concerning the Poisson's ratio, our calculated u=0.102 
for the pure BeH2 is much smaller than the experimental 
measurement on the amorphous phase. This is also re- 
ally understandable. In fact, it is well known that for the 
common materials that have much smaller shear moduli 
compared with the bulk moduli, their Poisson's ratio is 
close to 1/3. If the shear modulus, on the other hand, is 
much larger than the bulk modulus, materials will have 
negative Poisson's ratioi^. In the case of BeH2, our cal- 
culated results in Table I show that the shear modulus 
is comparable to the bulk modulus. Then, it is antici- 
pated that the Poisson's ratio of the orthorhombic BeH2 
has a small value, which is consistent with our calculated 
result of t'=0.102. The experimentally measured larger 
Poisson's ratic^ is therefore due to the use of amorphous 
sample. 

B. Electronic structure and charge distribution 

The calculated band structure of BeH2 is shown in 
Fig. 2. Our theoretical valence-band width is 6.5 eV 
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FIG. 2: GGA band structure of BeH2. The Fermi energy is 
set at zero. 



and band gap 5.5 eV, which is consistent with previous 
calculatioiiSi. As an illustrating comparison, in isoelec- 
tronic MgH2 the calculated valence-band width is 7.5 eV 
and band gap 4.2 eV-i2,, which implicitly suggests that 
the Mg— H bond is more ionic than the Be— H bond. Also 
one can see from Fig. 2 that the valence band maximum 
(VBM) and conduction band minimum (CBM) both ap- 
pear at X point in the Brillouin zone. By calculating 
the orbital-decomposed band charges inside the mufhn- 
tin spheres of BeH2 at VBM and CBM at X point, we 
find that the VBM has predominant H Is state charac- 
ter mixed with significant Be 2p contribution, while the 
CBM has obvious Be 2s and 2p features mixed with a 
little H Is contribution. 

The total electronic density of states (DOS) per unit 
cell for the orthorhombic BeH2 is shown in Fig. 3(a). For 
more clear illustration, the orbital-resolved partial DOS 
per atom in the unit cell are also shown in Fig. 3, see 
Fig. 3(b) for HI and H2 atoms, and Fig. 3(c) for Bel and 
Be2 atoms. From Fig. 3 the following prominent features 
can be seen: (i) H Is state hybridizes with Be 2s and 2p 
states in the whole range of the valence band. This is dif- 
ferent from alkali-metal hydrides. Due to the combined 
fact that there has only one valence electron with low 
ionization energy and the s— » p promotion energy is rela- 
tively large in alkali-metal elements, usually alkali-metal 
hydrides have entire ionic bonding in the sense that the 
valence electrons in alkali-metal atoms totally transfer to 
the H atoms. This clear ionic picture in alkali-metal hy- 
drides becomes blurred in alkaline-earth-metal hydrides. 
In the present case of BeH2 , in addition to a large ionic 
weight (see the following charge-distribution analysis), 
a distinct covalent component also exists in the Be— H 
bond, which is responsible for the stabilization of the low 
symmetric structure of Bell2. Prominently, this cova- 
lency is dominated by the hybridization of H Is and Be 
2p orbitals near the top of the valence band, while the hy- 
bridization of H Is and Be 2s orbitals is relatively weak, 
as shown in Fig. 3. Thus, the ionicity in the Be— H bond 
is mainly featured by charge transfer from Be 2s to H 
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FIG. 3: (Color online) Total and orbital-resolved local densi- 
ties of states for ground state BeH2. The Fermi energy is set 
at zero. 



Is atomic orbitals. Interestingly, in beryllium oxide the 
Be— O bond also exhibits strong covalency. As a result, 
bulk BeO has wurtzite structure instead of the usual form 
of rocksalt structure taken by other alkaline-earth-metal 
oxides; (ii) In the range of conduction band the DOS is 
mainly featured by Be 2s and 2p states, mixed with a 
little H Is states; (iii) Although there are three Be— H 
bond configurations (Bel-Hl, B1-H2, and Be2-H2), 
the DOS propertie of these chemical bonds are almost 
the same and there is no remarkable difference between 
Bel and Be2 (or between HI and H2) atoms. 

In order to gain more insight into the nature of Be— H 
chemical bonding in the orthorhombic BeH2 , we have in- 
vestigated the valence charge density distribution. Two 
contour plots of the valence charge density distributions 
are shown in Fig. 4. Here the contour plane in Fig. 4(a) 
is the plane composed by the three labeled atoms (1 to 
3) shown in Fig. 1, while the contour plane in Fig. 4(b) 
is the (001) plane (perpendicular to the c axis). It is 
these two planes that characterize the bond structure of 
BeH2. In both planes, hydrogen charge density shape 
is deformed toward the direction to the nearest-neighbor 
Be atoms. The charge density around H atoms is higher 
than that around Be atoms. This further indicates afore- 
mentioned significant ionic-type charge transfer from Be 
2s to H Is state. On the other hand, the charge densi- 
ties at the midpoint of Be— H bond is 0.46 e/A'^, which 
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FIG. 4: Valence charge density maps of BeH2 where (a) is 
the plane established by the three labeled atoms (1 to 3) as 
shown in Fig. 1 and (b) is the (001) plane. Here • stands for 
Bel, ▲ for Be2, o for HI and A for H2. The contour lines are 
drawn from 0.0 to 1.0 at 0.05 ev/A^ intervals. 



MgH2. In this way, H ionic radius is defined as the half 
of average nearest-neighbor H-H distance. Subtracting 
the H ionic radius from the average nearest-neighbor 
Bc-H distance gives the ionic radius of Be. Thus the 
ionic radius of H in our system is 1.131 A and 0.237 A 
for Be ions. After integration, we find 0.005 electrons 
around Be and 1.631 electrons around H. As a result, 
the valency of Be and H are represented as Be^-^^+ and 
For MgH2 our calculations give the valency of 
Mg and H as Mgi-^^+HO ™-. Clearly, our calculated 
valency shows that compared to the MgH2 system, 
charge transfer process is more prominent in BcH2 
system. Therefore, we would like to derive that the the 
definition of ionic radius by Yu and Lam still works for 
more covalent BeH2 hydride. 



is prominently higher than that in typical ionic crystals. 
This values is also much larger than that in isoelectroic 
MgH2 system. In fact, our calculation on Mgll2 gave 
charge density of 0.17 e/A^ at the midpoint of the Mg— H 
bond, which is comparable well with recent experimental 
datai^ of 0.2 e/A^ Thus again, the charge distribution 
reveals that hydrogen is more strongly bonded to Be in 
Bcll2 than to Mg in Mgll2. As a result, the calculated 
Be— H bond length (^1.37 A, see Table I) is remarkably 
smaller than the Mg— H bond length (1.93 A). Although 
hydrogen atoms do not form bond, for illustration, here 
we also present the charge densities at the midpoint in 
the line that links the nearest-neighbor H atoms, which 
are ~0.2 e/A^ for BeH2 and -0.15 e/A^ for MgH2. In- 
terestingly, the nearest-neighbor H-H distance in BcH2 is 
also larger than that in MgH2 by ~0.25 A. 

From above results, now it becomes clear that the 
bonding nature of hydrogen in BeH2, including the 
ionicity and covalency, is even more complex than 
in MgH2. Specially, because of complicated charge 
transfer process, which involves not only ionic but also 
covalent charge redistribution, one may wonder how to 
characterize the valency of H and Be in BcH2. Without 
knowing more appropriate way to valency, here we define 
and determine the so-called ionic radii of H and Be by 
using the method raised by Yu and ham^ in studying 



IV. CONCLUSION 

In summary, we have investigated the structural, me- 
chanical, electronic, and chemical bonding properties 
of BeH2 through first-principles DFT-GGA calculations. 
The calculated structural parameters are in good agree- 
ment with available experimental results. We have cal- 
culated the nine independent elastic constants and our 
results well obey the mechanical stability criteria. In par- 
ticular, the calculated C12 and C23 are negligibly small, 
and C13 is much smaller than Cn or C33. The bulk mod- 
ulus B, shear modulus G, Young's modulus E and Pois- 
son's ratio v have been obtained based on the knowledge 
of the elastic constants. The bonding nature of hydrogen 
in BeH2 have been fully analyzed in terms of the band 
structure, density of state, and valence charge distribu- 
tion. It has been shown that the Be— H bond displays 
a mixed ionic/covalent character. Here the ionicity is 
mainly featured by charge transfer from Be 2s to H Is 
states, while the covalency is manifested by prominent 
hybridization of H Is and Be 2p states. The elemental 
valency in BeH2 has been determined to be Be^-^^+ and 
By a systematic comparison, it has been shown 
that the chemical bonding nature of hydrogen in BeH2 is 
even more complex than in MgH2. 
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